Stability of Binary Mixtures in Electric Field Gradients 
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We consider the influence of electric field gradients on the phase behavior of nonpolar binary 
mixtures. Small fields give rise to smooth composition profiles, whereas large enough fields lead 
to a phase-separation transition. The critical field for demixing as well as the equilibrium phase 
separation interface are given as a function of the various system parameters. We show how the 
phase diagram in the temperature-composition plane is affected by electric fields, assuming a linear 
or nonlinear constitutive relations for the dielectric constant. Finally, we discuss the unusual case 
where the interface appears far from any bounding surface. 



I. INTRODUCTION 



The effect of gravitational and magnetic fields on the 
phase diagram of hquid mixtures is quite small in gen- 
eral due to the weak coupling of the field with the mix- 
ture's composition. The influence of electric fields was 
studied extensively, extensively in geometries where the 
field is uniform. Theoretical work predicted that for a 
binary mixture the upper critical solution temperature 
Tc is shifted upward by a small amount, of the order of 
milikelvinaii^. Experiments in low molecular weight liq- 
uids predominantly showed an opposite shift of the same 
magnitude^'^'^'^. Two exceptions are Reich and Gordon— 
who measured the cloud point temperature of a polymer 
mixture and Gabor and Szalai— who predicted a down- 
ward shift in the critical pressure of dipolar fluid mixtures 
in uniform electric fields. 

In a uniform field the shift in Tc is proportional to 
the square of electric field E'^ and (Pe/d(j)^ - the second 
derivative of the dielectric constant with respect to the 
mixture compositiona l "'^"'^^i^^'^'^i-'^'^'^^ . However, in realis- 
tic systems, such as colloidal suspensions and microflu- 
idic devices, the electric field varies in space due to the 
complex geometry. In such systems, electric fields on the 
order of lO^V/m naturally occur due to the small length- 
scales involved. Recently, we have shown that a homo- 
geneous mixture confined by curved charged surfaces un- 
dergoes a phase-separation transition and two distinct 
domains of high and low compositions appear— lil. The 
transition occurs when the surface charge (or voltage) 
exceeds a critical value. In this paper we examine in 
detail numerically and analytically the location of the 
transition. We also show how this transition affects the 
temperature-composition phase diagram. We find that 
near Tc the spatial variation of the field leads to a non- 
trivial modification of the stability lines. Our analysis 
shows that in nonuniform electric fields even a linear 
constitutive relation can lead to a substantial change of 
the transition (binodal) temperature. Lastly, we demon- 
strate that the phase separation interface can appear far 
from any of the surfaces bounding the mixture. 



II. THEORY 

Consider an A/B binary mixture confined by charged 
conducting surfaces giving rise to an electric field. The 
free energy of the mixture is 



F = I [/„ 



,r) + /e.(0,E)]dr, 



(1) 



where is the bistable mixing free energy density, given 
in terms of the A-component composition (0 < (/) < 1) 
and the temperature T. 

fes is the electrostatic free energy density due to the 
electric field E, given by 



(2) 



4> is the electrostatic potential. The negative sign cor- 
responds to cases where the potential is prescribed on 
the bounding surfaces. When the charge is given on the 
bounding surfaces the sign should be replaced by a pos- 
itive one^ii^i^. e is the mixture's permittivity, and is a 
function of the composition (j>. Note that in order to iso- 
late the electric field effect we do not include any direct 
short- or long-range interactions between the liquid and 
the confining surfaces. The equilibrium composition pro- 
file (/)(r) and electrostatic potential ^p{r) are given by the 
extremization of the free energy with respect to (f) and 
tfj22i^. The resulting Euler-Lagrange equations are 



SF 
6£ 



Sfrn 1 cie(0) 



(V^)2 - M = 0, 



V • (£(0)Vi/;) = 0. 



(3) 
(4) 



Eq. ([4]) is simply Gauss's law. Notice that the relation 
£(0) couples these nonlinear equations. In the canonical 
ensemble is a Lagrange multiplier adjusted to satisfy 
mass conservation: 



J - 0o] 



dr = 0, 



(5) 



where 4>o is the average composition. In the case of a sys- 
tem in contact with a matter reservoir (grand-canonical 
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FIG. 1: The three model systems, (a) A single charged spher- 
ical colloid with surface charge density a and radius (b) 
A charged wire with surface charge density a and radius Ri, 
or two concentric cylinders with radii -Ri and 7?2. (c) A wedge 
comprised of two flat electrodes with an opening angle f3 and 
potential difference V. Ri and R2 are the minimal and maxi- 
mal values of the distance r from the imaginary meeting point 
of the electrodes. 



ensemble), the chemical potential is set by the reservoir, 
fi = iJ-{4>o) where 0o is the reservoir composition. 

In order to simphfy the solution of Eqs. Q-®, we 
consider the three simple model systems shown schemat- 
ically in Fig. [T] The first one is a charged isolated 
spherical colloid of radius Ri and surface charge den- 
sity cr, immersed in an infinite mixture bath. In this 
system, spherical symmetry dictates that (f) = 4>{r) and 
E = E(r) where r is the distance from the colloid's center. 
Since the colloid has a prescribed charge we can integrate 
Gauss's law and obtain an explicit expression for the elec- 
tric field: E(r) = aRl/{e{(j))r'^)r. The second geometry is 
a charged wire of radius i?i and surface charge density cr, 
coupled to a reservoir at r — > 00. Alternatively, we may 
consider a closed condenser made up of two concentric 
cylinders of radii i?i and i?2- In both cases, we readily 
obtain the electric field E(r) — aRi/(e{4>)r)r, where r is 
the distance from the inner cylinder's center. The last 
system is the wedge condenser, made up from two flat 
electrodes with an opening angle /3 and a potential dif- 
ference V across them. Solution of the Laplace equation 
gives E(r) — {V/ (3r)6, where r is the distance from the 
imaginary meeting point of the electrodes and is the az- 
imuthal angle. In this geometry (p = (f){r) and therefore 
E • Ve = 0. The explicit expressions for the electric field 
in all three systems outlined above decouple equations 
© and Q. 

We will show that Eq. ([3]) leads, under certain con- 
ditions, to a phase-separation transition. This transi- 
tion is independent of the exact form of fm, and can 
be realized as long as fm is bistable and the dielectric 
constant e depends on the composition </>. In order to 
be specific, we will consider the mixing free energy de- 
rived from the Flory-Huggins lattice theory, with lattice 
site volume vq- We consider the simple symmetric case 
where each component occupies N successive lattice cells. 
Simple liquids have = 1, while polymers have A'^ > 1 
monomers. The mixing free energy density is then given 
by fm = ksTfm/Nvo, where 

fm = plogW + (1 - 0) log(l - 0) + iVx0(l - 0). (6) 

is the Boltzmann constant, and x ~ 1/T is the Flory 
interaction parameter™. We limit ourselves to the case 



where x > Oj leading to an Upper Critical Solution Tem- 
perature type phase diagram in the (f> — T plane. In 
the absence of electric field, the mixture is homogeneous 
above the binodal curve (f>t (T) , and phase separates into 
two phases having the binodal compositions 0t below 
it. Below the binodal curve, but above the spinodal, 
given by (j)s{T) = (1/2) [1 ± ^1 - 2/{Nx)], the mixture 
is metastable. The binodal and spinodal curves meet at 
the critical point {(j)c, {Nx)c) = (1/2,2). The transition 
(binodal) temperature Tt for a given composition is given 
by Tti4.) = {Nx)cTc [log(0/(l - </.))/(20 - 1)]-^^. 

Using the expressions given above for the electric field, 
we write the generalized composition equation valid for 
cylindrical and spherical geometries: 



/;„(0)-^xM.ef|^f-"-A = o. 



(7) 



Here, f = r/Ri is the scaled distance and e = s/sq, with 
Eq the vacuum permittivity. Where 



(8) 



is the dimensionless field, and n is an exponent charac- 
terizing the decay of £^^: n = 2 for concentric cylinders, 
and n = 4 for spherical colloid. For the wedge geometry 
we find: 



di 



where 



Afi^ksT^Rl ■ 



(9) 



(10) 



and n — 2. Mgc and are dimensionless quantities 
measuring the magnitude of the maximal electrostatic 
energy stored in a molecular volume compared to the 
thermal energy. The second term in Eqs. ([7]) and ^ is 
the variation of the electrostatic free energy with respect 
to 0, and is only present when the mixture components 
have different permittivities. 

The constitutive relation e{4>) is a smooth function of 
(j). Experiments show that the curve can be slightly con- 
vex or concave, and is dominantly linea r They 
are mostly in agreement with Clausius-Mossotti and 
Onsager-based theories for the dielectric constant^^. 
Thus, for a mixture of liquids A and B with dielectric 
constants ia and £{,, respectively, the experiments yield 
a polynomial relation in the form 



e{(f>) ~ £6 -I- -I- £"0^ -t- 



(11) 



We start by focusing on a linear relation, namely e" — 0; 
in this case e' — Ae = Sa — st- Even in such a simple case 
it turns out that a phase separation transition occurs, 
in contrast to the Landau mechanism which relies on a 
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FIG. 2: Graphical solution of Eq. (|12p for an open wedge at 
T < Tc and a symmetric mixture. Solid curve is f'm{4')- Its 
roots are the transition (binodal) compositions. The intersec- 
tion between ,f'm{4') a-nd the horizontal dash-dotted line ji{f) 
gives the composition (t){f). If p.{Ri) is at p.a, the profile (/)(f) 
varies smoothly, but if p.{Ri) — fib, has a discontinuity. 



FIG. 3: The three types of equilibrium profiles (f)(f) for a sys- 
tem of two concentric cylinders. Dash-dot line: T = l.OOSTc 
and Msc = 0.04: above Tc the profile is smoothly varying. 
Dashed line: T = 0.995Tc, M^c = 0.01 smaller than the crit- 
ical value for demixing. Solid line: same T, but Msc = 0.04 
is large enough to induce phase separation marked by an in- 
terface at f = 7?. We took an average composition (j>o — 0.41. 
Here and in other figures, Ri = 1, R2 = 5, Sa ~ 5, and it ~ 3. 



nonvanishing e" . After investigating linear relations we 
examine how our results change when e((/i) has a positive 
or negative curvature, by allowing for e" ^ 0. In this 
case e' is different from Ae. Higher order terms in the 
expansion Eq. (|lip are not expected to change the results 
qualitatively, since they do not affect much the curvature 
of e(0)^^^. 



III. RESULTS AND DISCUSSION 

Before we present the numerical solutions of Eqs. ([7]) 
and ([H]), it is illustrative to consider a graphical solution 
for a wedge condenser. Recall that in the absence of field 
it is assumed that T is above the binodal temperature. 
We rewrite Eq. ^ as: 



/^W = S5 



(12) 



Mo 



where /xq is the dimensionless reservoir chemical poten- 
tial. At a given temperature, the intersection of f^i{(j>) 
and the horizontal line /i(r) gives the local composition 
4>{f) (Fig. [2]). When f — > 00, fl{r) — > p,o and the compo- 
sition is = 00 1 corresponding to a homogeneous phase. 
For simplicity we consider < (j)c- As f decreases, /2(f) 
(and hence (/'(f)) increase until they attain their maximal 
value at Ri. Above Tc, the free energy is always convex, 
is a monotonic function of (p, and the composition 
profile (/'(f) is hence continuous. However, below Tc, fm 
is bistable and is sigmoidal. In this case there are 



two possible scenarios shown in Fig. [5] If, for example, 
/2(i?i) = jla, the composition profile varies smoothly. If, 
on the other hand, /i(i?i) — jib, there is a discontinuity 
in the profile since there is a radius R > Ri where the 
value of (j) can "jump" from high to low values. Below 
Tc, there is a range of radii, or compositions, where the 
discontinuity in ^(f) can occur. The equilibrium profile 
(j){f; R) is the one that minimizes the total free energy 
integral F = J f[(j){f;R)]dr. These conclusions also hold 
for Eq. 0. 

In Fig. [5] we show the composition (j)s defined by the 
relation f^{(t)s ) = finis's)- Clearly, this is the minimal 
composition for which exist more than one solution to 
Eq. The role of this special composition will be 

discussed later. 

The three typical composition profiles are shown in 
Fig. El Above Tc (dash-dot line), (f>{r) varies smoothly 
due to the dielectrophoretic force, whereby the high-e liq- 
uid is drawn into the strong electric field region. Below 
Tc (dashed line), at a temperature where the field- free 
mixture is homogeneous, if M is small the profile 4'{f) 
is again smoothly decaying, exhibiting the same dielec- 
trophoretic behavior. However, if at the same temper- 
ature, M is increased, either by adding charge to the 
surface or by increasing the curvature (smaller Ri), we 
arrive at a critical value, denoted M*. Above it, a phase- 
separation transition occurs. This is shown in the solid 
line of Fig. [31 where the mixture consists of two coexist- 
ing domains separated by an interface at f = i?. 

The typical value of charge/ voltage required for demix- 
ing can be estimated from the value of M being in the 
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FIG. 4: Location of the demixing interface J? as a function 
of M, (a) for three average compositions 0o in a closed cylin- 
drical system (lines with symbols) and for 0o = 0.28 (dashed 
line) and (po — 0.3 (dash-dot line) in an open system. Curves 
do not coincide at i? = 1. (b) The same, in the closed "wedge" 
geometry. We took T = 0.96Tc. 



range M ^ 0.001 — 0.1-12.. Consider a colloid of radius 
Ri ^ l/im placed in a mixture having a molecular vol- 
ume Nvo ^5x 10-"m3 and ~ 300K. Then, the typ- 
ical demixing charge is of the order of 10^ — lO'^e charges 
(surface voltage is 1 — 1001^). It scales linearly with 1/iV: 
in a polymer mixture the confining surfaces require N 
times smaller charge compared to molecular liquids in 
order to induce phase separation. 



A. The phase-separation interface 

At the critical value of AI, a sharp interface first ap- 
pears separating coexisting regions of high- and low-0 
value. If the average composition (jjQ is smaller than (j^c, 
the interface appears at r = When we further in- 
crease M and supply more electrostatic energy to the 
system, dielectrophoresis leads to an increase in the size 



of the high composition domain. Thus, the location of 
the separation interface R increases. Fig. 2] shows how 
R varies with M at a constant temperature in the con- 
centric cylinders and wedge systems. Notice that as (po 
approaches the binodal composition (0t ~ 0.33), R is 
larger at the same M . It also grows more rapidly with 
increasing (j)o. Indeed, when the binodal is approached, 
the mixing free energy barrier is smaller. Secondly, R 
grows faster in an open system than in a closed one. 
This is because in a closed system the energy penalty 
in fm grows faster than the energy gain in fes- Material 
conservation gives the maximum value of R, Roo, given 

by 



Ri 



(13) 



This M ^ oo limit is physically unattainable and is pre- 
empted by dielectric breakdown. 

Note that the typical values of M in the cylindrical and 
spherical cases are an order of magnitude larger than in 
the wedge condenser, see Fig. |4] (a) and (b). Indeed, 
in the spherical and cylindrical symmetries, E is paral- 
lel to V0: the dielectrophoretic force, (proportional Ae), 
has to be large enough to overcome the energy penalty 
associated with dielectric interfaces parallel to E (pro- 
portional to (Ae)^)i^. On the other hand, in the wedge 
condenser E is perpendicular to V^, and the required 
dielectrophoretic force for demixing is correspondingly 
smaller, leading to smaller values of M^. This could be 
seen by comparing the electrostatic terms in Eqs. ^ 
and differing by a factor of e((/))~^ ~ 0.1, which Msc 
has to compensate for in order for the values of fes to be 
equal. 

Alternatively, an increase in T at constant electric field 
decreases xM ~ M/T and decreases R. Fig. [5] shows 
how an increase in T shrinks the high-(/) domain and de- 
creases i? in a closed and open cylindrical system. In 
Fig. [5l when R — \ the temperature is that for which 
M = M*. The effect of temperature is much more pro- 
nounced in an open system: R tends to infinity when ap- 
proaching the binodal temperature (not shown). On the 
other hand, R is finite when approaching the binodal in a 
closed system. Its maximal value is larger when l^o — (/"d 
is smaller, because then the mixing free energy difference 
between low- and high-c/i values is reduced. In Fig. [5] (a), 
R appears to be linear simply because T changes over a 
small interval. 

One can estimate the value of 4> at the demixing in- 
terface in an open system. At the interface there is a 
"jump" in (/)(f) from (j>i to (/)2 > (pi. Let us denote by 
4>{R) the upper interface composition - 0(i?) = (t>2- The 
conditions for a "jump" in the wedge geometry are: 



AeM^ 



Mo = 0, (14) 



/™(02) - 



i?2 

iVx^-Ao^O, (15) 
,{(t)2)-Nx . (16) 



5 




0.935 



0.945 0.955 

T/T, 



0.965 




0.94 



0.95 0.96 0.97 

TIT, 



FIG. 5: Location of the demixing interface _R as a function 
of T /Tc, (a) for three average compositions in a closed cylin- 
drical system: 0o = 0.28 (squares), 0.3 (circles) and 0.32 
(crosses) with Msc = 0.04. When R — Ri, the tempera- 
ture corresponds to Msc = M^^ = 0.04. R grows as T/Tc is 
reduced until it attains it maximal value at the binodal tem- 
perature (dashed line for each value of 00 ) (b) Same, in an 
open cylindrical system with Msc = 0.004. Here, _R ^ cxd at 
the binodal. 



The first two equations define the local solutions of Eq. 
(fT2)) . and the third one is the condition that a high com- 
position is favorable: f{4>2) < f{4'i)- The true value 
of R is the one that gives the global minimum of the 
free energy. Putting Eq. I|14p in Eq. and using 

Ao = /m(</'o), we get 



fm{4>2) - /m(0l) 



(17) 



The right-hand side of Eq. pT)) is maximal when the 
transition occurs from 0i 
We therefore denote A/^. 
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(18) 
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FIG. 6: (a) The "differentiating" curve - (j)' (dashed line) for 
an open wedge. Above it, when phase separation occurs, the 
composition at the interface jumps from (/)s^' to cf>s- Below (f)^ 
the front composition jumps from ^i^' < cf>\ < <j>t to (j>2 > 4's- 
Examples of this behavior are shown in (b) and (c), where 
composition profiles for points above and below (j>^ with two 
values of A/ are given, showing this behavior is independent 
of M. The inset in (a) is a blowup showing the location 
of (b) and (c) in the <j)o — T plane. In (b) <^o = 0.38 with 
Mw = 4 X 10"" and Mw = 6 x 10"", and in (c) = 0.39 
with Mw = 1 X lO"'' and M„ = 1.5 x 10"^ In (b) and (c) we 
took r = 0.987rc. 



If the inequality 



/m(0l) - /m 



> A/, 



(19) 



(21 

holds, the transition must be at (j)i — (pi , since for larger 
values of 0i the right-hand side of Eq. (fT7|) is smaller 
while the left-hand side is larger, so a higher composition 
is surely favored. The equality sign in Eq. p9p corre- 
sponds to the maximal average composition 0o for which 
this equation holds. 

The locus of such compositions is the "differentiating 
curve" - (t)\ shown in the dashed curve of Fig. [5] (a). 
When the zero- field point in the phase diagram {(l>o,T) 
IS above the ^t(T) curve, Eq. (fT9)l holds, and the com- 
position at the interface jumps from c/ii^'' to (pi^^. When 
{(pQjT) is below (j)^{T) the upper interface composition 
<j>{R) is between 0s and (j>f Eq- shows this result is 
independent of M. An example of this situation is given 
in Fig. [6] (b), where the compositions 4>i and 02 are the 
same for two values of M. Fig. [5] (c) shows the inter- 
face compositions are independent of AI also when 0o is 
larger than cj)^ . 
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FIG. 7: Composition at the demixing interface i? as a function 
of M for a closed wedge (crosses) and concentric cylinders 
(circles). The dashed line is the binodal composition. Here 
(l>o = 0.3 and T = 0.96Tc. 



When the discontinuity in the composition profile oc- 
curs at 02 = (/"s, we can invert Eq. © to get R oc M-^/^ 
for the open wedge system. In particular, above the dif- 
ferentiating curve we get in the Flory-Huggins model 



R = 



fMs) - Mo 



1/2 



(20) 



Thus, R varies linearly with the wedge potential V . 

In the other geometries /g^ depends on 0, and this 
influences the value of pairs 0i and 02- However, a very 
good approximation, valid when Mgc is not too large and 
T is not too close to Tc, is that the transition remains 
from (pi — (ps to 02 — 4>s ■ One can then repeat a similar 
derivation and obtain 



f' 

•I m 



A/. 



SC. max 



,max; 



/,(2) 



a(2) 



(21) 
(22) 



Using these equations one can determine the differenti- 
ating curve for cylindrical and spherical geometries. 

The situation is different in a closed system, as Fig. 
[7] shows. In the wedge, the demixing interface occurs at 
the binodal composition, 0t, irrespective of M. In the 
cylindrical and spherical systems, 0(i?) is lower than but 
close to (f>t, and decreases when M grows. The qualita- 
tive explanation is as follows. In the wedge geometry /^^ 
only adds a constant to /', and the binodal compositions 
remain the only pair of solutions of Eq. ([9]) that have the 
same mixing free energy /,„. Hence, the mixing free en- 
ergy penalty is minimized when the transition is at the 
binodal compositions^^. This also explains why in the 
wedge geometry 4){R) is independent of M . In the cylin- 
drical and spherical geometries on the other hand, Jig 
affects (t){R), resulting in a value of (j>{R) smaller than 





FIG. 8: (a) Stability diagram of a spherical colloid with Msc = 
0.04 (thick solid line) and Msc = 0.08 (dashed line). Dash-dot 
curve is 4>^ [see Fig. |6ja)] and thin solid line is the binodal 
Tt{(j)o) ■ (b) Stability diagram showing <j)*{T) for concentric 
cylinders with Msc = 0.04 (thick solid line) and Mac = 0.08 
(dashed line). 



<f)t- This reflects the fact that dielectrophoresis favors 
high values of 0. Since fes oc M, larger values of M lead 
to lower values of (f>{R). 



B. Stability diagrams 

One can also set M constant and for a given electric 
field draw the a stability curve (j)*{T) in the (j)o — T plane, 
see Fig. m (j)*{T) is defined such that below it phase 
separation occurs, while above it composition profiles are 
smooth. Fig. [5] (a) and (b) show the stability diagram of 
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FIG. 9: Composition profiles above the kink temperature for 
an isolated spherical colloid. The profiles change from smooth 
(dash-dot and dashed lines) to discontinuous (solid line) with 
a discontinuity at finite value of 7?: Ri < R < R2, when 0o 
increases at constant T and AI (see inset). 



a spherical colloid and concentric cylinders, respectively. 
In both diagrams, an increase in M increases the unstable 
region. For the same value of M, in a closed system 
the phase separation region is smaller than in an open 
one, because the mixing energy penalty makes it more 
difficult to induce phase separation. The range of values 
of (po that are unstable in nonuniform electric fields grows 
when T increases (but still T < Tc). For low values 
of T, there is a significant difference between open and 
closed systems: in open systems, if M is large enough 
the stability curve tends to cj)^ [see Fig. EJa)], whereas 
for closed systems the stability curve tends to (j)t where 
demixing is spontaneous. 

In both parts of Fig. [8l there is a kink in all the curves 
0*(T) at a temperature we denote Tk^i- In the Flory- 
Huggins model and for spherical colloids and concentric 
cylinders, the second derivative of the free energy is 



d^fn 



2NxM, 



jde/dct))^ 



(23) 



The second term in this equation is the positive electro- 
static contribution f'J^. It is clear that when an electric 
field is present, even at T < Tc, the electrostatic contri- 
bution can lead to a positive value of f"{4>, f) and phase 
separation cannot occur. 

When T < Tk^i, f" is negative for all values of f, and 
the phase separation interface appears first at i? = 
However, when T > T^ i, for a given value of 0o, a special 
radius Rc{T, (j)Q) exists. This is the largest value of f for 
which /" in Eq. I|23p can be negative. In this case, the 
demixing interface appears first at R = Re- An example 
of this behavior is shown in Fig. [HI at constant Tc > 
T > Tfe^i, at points A and B (0o < 0*) 0(^) is smooth, 
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FIG. 10: Surface composition (^(i?i) when approaching the 
binodal, (a) for three temperatures in a closed cylindrical sys- 
tem, (b) Same, for an isolated charged cylinder. We took 
M,c = 0.04. 



similar to (j){r) above Tc- However, at point C {(f>o > 4>*) 
(j){r) has a discontinuity at i? > Re- As the critical point 
is approached, Rc ^ 00 and (jf approaches the critical 
composition. 

The kink temperature Tfc i is given by setting Rc = i?i 
and can be obtained from solution of: 



(24) 



We stress that this result is independent of the exact 
form of mixing free energy. Notice that for a wedge, the 
electric field has no effect on the convexity of the free 
energy, and one finds Tk^i — Tc (see Fig. [TTt . 

The surface composition when approaching the 

binodal at constant T and M is given in Fig. [TOl When 
T < Tfc i the surface composition has a discontinuity at 
(j>o — (p* (dashed and dash-dot lines) and the value of 
<j>{Ri) becomes larger than (pc- When T > Tk.i the sur- 
face composition varies smoothly (solid lines). The dis- 
continuity in (t>{Ri) occurs at lower values of 0o in open 
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C. Quadratic constitutive relation 
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We now examine how the stabihty diagram changes 
if the dielectric constant has a quadratic dependence on 
composition: e" 7^ in Eq. (dJ). For simpHcity, we 
treat the wedge system where a hnear constitutive rela- 
tion means that T^^i = Tc- In the Flory-Huggins model, 
the conditions in Eq. l(24|) give 



Tk.i 



(25) 



where we used R\ — \. 

Note that e" can be positive or negative. When e" < 0, 
/"^ is positive and we return to the same behavior we saw 
for spherical and cylindrical systems with linear relation 
e(0). On the other hand, if i" > then /"^ is negative 
and phase separation is possible above Tc- The stability 
diagram of a wedge with a quadratic constitutive rela- 
tion is shown in Fig. [TTJ In this figure, arrows labeled 
A-E indicate the variation of 4>q at constant T in differ- 
ent areas of the stability diagrams. For each arrow the 
location for which the interface first appears is given in 
the caption of Fig. [TlJ For the data in Fig. [TTl the 
kink temperature is given by Tk/Tc = 1 ± 0.0064, de- 
pending on the sign of e" . Taking Tc w 300K, we have 
AT{(t)*{Tk)) = Tk-Tc ~ ±2K. This change in is 
two orders of magnitude larger than the corresponding 
change in uniform electric fields. Note that in most of 
the phase space the displacement of the transition tem- 
perature due to a nonvanishing e" is much smaller than 
that due to e'; in spatially nonuniform fields far from 
the critical composition, the dcmixing transition is well 
described by a linear constitutive relation 



D. Demixing transitions for (f) > 



FIG. 11: (a) Stability diagram for an open wedge with 
Mw = 4 X 10~^ and three constitutive relations Eq. (Illf) 
with £a = 5 and ib ~ 3. Thick solid line is the transition 
(binodal) curve Tt(c/>). When e" = (dash-dot line) there is 
no phase separation above Tc- For e" — 1.6 an unstable re- 
gion exists above Tc- In contrast, when i" — —1.6 is negative, 
the unstable region is suppressed. <j){f) changes from smooth 
to discontinuous at _Ri along lines A and B, and at Rc > Ri 
along C and D. (b) The same for a closed wedge system. Note 
that the stability curve extends to compositions larger than 
(j>c- 4>{r) changes from smooth to discontinuous at _Ri along 
A and H, a.t Ri < Rc < R2 along C and D, and at R2 along 
E. Here Tk,i/Tc = 1.0064. 



systems compared to closed ones; open systems are less 
stable than closed systems. When T increases at a given 
value of (j)o, the surface composition decreases since mix- 
ing is favored at high temperatures. 



Since the electric field breaks the symmetry of the free 
energy with respect to composition {(f> 1 — 0) , the full 
stability diagram is asymmetric with respect to cj) — ipc- 
Fig. [m (a) shows that in an open system, phase separa- 
tion does not occur when (j)o > ipc- In Fig. [TT] (b) there 
are unstable compositions (jjo such that 4>q > (pc- This is 
an important feature of the stability diagram in closed 
systems. Here, the dielectrophoretic force that pulls the 
high dielectric constant liquid toward the electrode cre- 
ates a depletion in the region where electric field is low, 
and phase separation will occur if the composition at this 
region is close to the binodal composition. The stability 
curves in Fig. [TT] (b) show that higher values of potential 
or charge are required for demixing when 0o > (pc- In- 
deed, the ratio of electrostatic energies at the inner and 
outer radii is ~ R^/Rf- 

Examples of this behavior are shown in Fig. [121 
When (j) > (j>c, and M is small (dashed line) (/)(f) varies 
smoothly, exhibiting the usual dielectrophoretic behav- 
ior. However, if M is sufficiently large, that is M > M*, 
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0.65 



0.55 



0.45 




FIG. 12: Equilibrium profiles (/!>(f) for a closed wedge sys- 
tem with an average composition <j)o = 0.58 larger than (pc- 
Dashed line: M„ = 1 x 10~^ smaller than M^, and (f>{r) is 
smoothly varying. Solid line: M„ = 4 x 10~^ is large enough 
to induce phase separation. Dash-dot line: Mw = 8 x 10~^, 
the separation interface moves to a smaller radius. We used 
r = 0.993Tc. 



the long range dielectrophoretic force gives rise to a 
phase-separation transition near i?2 (solid line). If M 
is further increased the phase-separation interface moves 
to smaller radii (dash-dot line). 

Fig. [T3] shows how R varies with M and T in the wedge 
geometry for (j) > (j)c- At M*, the interface appears where 
the field is minimal, i.e at R2 = 5, see Fig. [13] (a). An 
increase in M decreases R and the volume of the high-(/) 
domain. Fig. [T3] (b) shows that at fixed M, R decreases 
with T and attains its minimal value at Tf . The minimal 
value of R is given by Eq. (fT^ . 

When e" ^ 0, the stability curves in Fig. [TT] (b) have 
another kink at > (/)c at a temperature we denote Tk^2- 
From the previous discussion, Tk^2 is obtained by replac- 
ing Ri by i?2 in Eq. (|24p . In the wedge geometry we 
thus find 



~t7 



= 1 



M^e" 
Rl 



(26) 



Hence, — Tk^2 is smaller than — T^^i by a factor of 



Rl 



IV. CONCLUSIONS 

We present a systematic study of electric field induced 
phase-separation transitions in binary mixtures on the 
mean field level. The behavior described by us comple- 
ments the findings of Onuki and co-worker o^'^°i^^ and 
Ben-Yaakov et al. who mainly focused on effects above 
the critical temperature or on miscible liquids. The tran- 
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» 0^4.9 




0.9915 0.9925 0.9935 0.9945 

T/T, 



FIG. 13: (a) Interface location R for a closed wedge system 
with an average composition (po = 0.58 > as a function of 
Mw at T = O.ggSTc. The minimal value of _R is = 3.86 
given by Eq. (b) _R as a function of T at M„ = 4 x 10"^. 

Dashed line is Tt. 



sitions should occur in any bistable system with a dielec- 
tric mismatch sufficiently close to the transition tempera- 
ture, e.g a vapor phase close to coexistence with its liquid. 
The differences between closed and open systems and be- 
tween the three geometries are highlighted. The stabil- 
ity diagrams in the temperature-composition plane are 
given. These diagrams show that the change in the tran- 
sition temperature is much larger in nonuniform fields 
than in uniform fields. We describe the special temper- 
atures Tfc^i and Tk^2 where the stability diagram has a 
"kink" . It should be emphasized that the location of the 
phase separation interface is not restricted to the vicinity 
of the confining surfaces, and can appear at a finite loca- 
tion R in the range of temperatures Tk^i < T < Tk.2 as 
described above in Fig. [TT] In other geometries, too, the 
demixing interface can be created far from any surface, 
e.g in quadrupolar electrode array. 

In order to test our predictions we suggest the follow- 
ing experiments. Consider a wedge capacitor partially 
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immersed in a binary liquid mixture near the coexis- 
tence temperature. Neutron reflectometry may be used 
to probe directly the composition profile^ and to deter- 
mine the critical voltage for demixing. Here we expect 
R to be proportional to the voltage difference across the 
electrodes. A different experiment can be realized by 
suspending a conducting wire in a dilute vapor phase of 
a pure substance near the coexistence temperature with 
the liquid. When a potential V is applied to the wire, 
its fundamental frequency of vibration perpendicular to 
its length changes from /o to fv because a liquid layer 
condenses around the wire. We expect that the ratio 
(/o//v)^ should be a linear function of V^. Alterna- 
tively, one can measure directly the change in the wire 
mass Am using a microbalance - Am should depend lin- 
early on V^. In these experiments, care must be taken to 
subtract wetting or confinement effects due to the elec- 
trodes. 

Capillary condensation has been described^Si^i^^'^i^ 
for colloidal suspensions in binary mixtures. According 
to our work, charged colloidal suspensions can floccu- 
late or be stabilized due to the formation of liquid layers 
around the colloids. The liquid layer may also influence 
the interaction between a colloid and a flat surface both 
far and close to the critical point^^. 

Our results may also be measurable in Surface Force 
Balance and Atomic Force Microscope experiments. In 
such apparatus, a bridging transition has been ob- 
served in binary mixtures^^ and explained by capillary 
forces^ii^. Since in most cases the surfaces are charged. 



a capillary bridge could be the result of the merging of 
two field-induced layers, leading to long range attractive 
forces between the surfaces. 

In the current work we neglect a (V(/))^ in the mixing 
free energy /m, since the system size was relatively large 
and the electrostatic energy acts throughout the whole 
system volume, and therefore is dominant We have 
verified that inclusion of such term leads to smoothing of 
composition profiles but otherwise to no other noticeable 
changes in the figures presented. In order to isolate the 
electric field effect, we have also not considered any di- 
rect short- or long-range interactions with the confining 
surfaces. In a future study it would be very interesting 
to relax these assumptions and to look at ever smaller 
systems. Here we are intrigued by the possibility to find 
qualitative, and not just quantitative, differences from 
our current profiles and diagrams. 
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